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Abstract 

We introduce the Umbral calculus into Clifford analysis starting from the abstract of the Heisen- 
berg commutation relation a;] = id. The Umbral Clifford analysis provides an effective framework 
in continuity and discreteness. 

In this paper we consider functions defined in a star-like domain O C R" with values in the 
Umbral Clifford algebra C'i'a,^ which arc Umbral polymonogciiic with respect to the (left) Umbral 
Dirac operator D' , i.e. they belong to the kernel of (Z)')*^. We prove that any polymonogenic function 
/ has a decomposition of the form 

f = fi+x'h + --- + {x')''-^fk, 

where x' = Xiei + - ■ • + x'„en and fj,j — l,...,k, are Umbral monogenic functions. As examples, this 
result recoveries the continuous version of the classical Almansi theorem for derivatives and establishes 
the discrete version of Almansi theorem for difference operator. The approach also provides a similar 
result in quantum field about Almansi decomposition related to Hamilton operators. 

Some concrete examples will presented for the discrete analog version of Almansi Decomposition 
and for the quantum harmonic oscillator. 

Keywords: Umbral Dirac operator, Almansi's decomposition theorem, Umbral polymonogenic func- 
tions, Quantum Harmonic Oscillator. 



2 



1 Introduction 



We introduce the Umbral Calculus into the ClifFord analysis. This provides a framework unifying conti- 
nuity and discreteness. We verifies this phenomena by establishing the Almansi decomposition in Umbral 
Clifford analysis. 

The essence behind continuity and discreteness is Umbral Calculus. Umbral Calculus has many appli- 
cations in combinatorics, special function theory, approximation theory, statistics, probability, topology 
and physics (see [1]). Umbral calculus ([H [3]) can be viewed as an abstract theory of the Heisenberg 
commutation relation [P, M] = id. In quantum mechanics P is the derivative and M the coordinate 
operator. The unification of the Umbral structure can be seen by taking the derivative P to be the usual 
derivative and the divided difference. 

In the Umbral Clifford analysis, we take the delta operator Ox^ as momentum operators and take x'j = 
^ (^XjO'xj + O'xjXj^ as position operators, where O'xjfix) — O^j {xjf{x)) — XjOxjf{x) is the Pincherle 
derivative. They satisfy the Heisenberg- Weyl relations 

[Ox, , OxJ = = [x'j,x'k], [Ox, , Xfe] = Sjkid. 
The Umbral Clifford analysis is a function theory for null solutions for the Umbral Dirac operator 

n 

D' ■.^J2^,Ox,. 

where ei, . . . ,e„ are the Clifford Algebra generators subject to the anti-commuting relations 

e^Gfe + ekBj = -26 jk- 

By taking Ox^ to be the partial derivative we recovery the standard Dirac operator D = ^j'^xj ■ 

The Clifford analysis is a refinement of harmonic analysis and a generalization of complex analysis to 
higher dimensions (see [H [5l [6l [7] ) • In the spherical harmonics, the Fischer decomposition of polynomials 
plays a central role. The Almansi decomposition generalized the Fischer decomposition from polynomials 
to polyharmonic functions. The polyhamonic theory has many applications in elasticity 9J, radar imaging 
[TU] , and approximation [TTl [T^] . 

Almansi's Theorem (cf. [13[ I14j ) /// is polyharmonic of degree k in a starlike domain with center 0, 
then there exist uniquely defined functions /i, • ■ • , fk, each harmonic in such that 

fix) = /i(.t) + \x\'f2{x) + ■■■ + |x|2(fc-i)/fe(a;). 



One can find important applications and generalizations of this result in the case of several complex 
variables in the monograph of Aronszajn, Creese and Lipkin, |14j, e.g. concerning functions holo- 
morphic in the neighborhood of the origin in C". More recent generalizations of Almansi's Theorem can 
be found in [miinillllllH]. 

In [19j . Malonek and Ren obtained the Almansi decomposition in Clo.n- In this paper, we establish the 
Almansi decomposition in the setting of Umbral Clifford analysis. As a corollary obtained a discrete 
version of Almansi decomposition. The approach also provides a similar result in quantum field about 
Almansi decomposition related to Hamilton operators. 
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2 Umbral calculus revised 



2.1 Shift-Invariant Operators 

Let X = {xi,X2, ■ ■ ■ ,Xn) e K". A monomial over x is the product x" = x"^X2^ ...a;"". Here and 
elsewhere a, (3, . . . denote indices of nonnegative integers. A polynomial is a finite linear combination of 
monomials. We denote the ring of polynomials over x by R[x], d^j '■= the partial derivative with 

respect to Xj and dx {dxi , , ■ ■ ■ , dx,^ ) the gradient operator. 
In the sequel, it will be helpful to introduce the shorthand notations 

• ^Xi ^X2 ■ ■ ■ Xn ' 

a! = q!i!q;2! . . . q;„!, 



a ) a!(/3-a)!' 



Set y ■ dx = X]'i=i Vj^xj and define the shift operator 



T, 



f{x) = exp(j/ • dx)f{x) = f{x + y). 



An operator Q is shift-invariant if it commutes with the shift operator Tj,, when acting on M[x], i.e. 

[Q,r^]p(x) = (1) 

holds for all p S ]R[^] and y £ M", where [a, b] = ab — ba denotes the commutator bracket between a 
and b. In the language of co- algebra, the shift operator Ty is nothing else than the co-product in a Hopf 
algebra over M.[x], but we shall not discuss this here, see e.g. [20] (Section V.3) for more details. 

Furthermore the following multivariate expansion formula holds |21| : 

Theorem 2.1 A linear operator Q : R[x\ — > R[x] is shift-invariant if and only if it can be expressed (as 
a convergent series) in the gradient dx, that is 

oo 

Q=Y1 

|a|=0 

where Oa — ^4- 

L J 1=0 



We will denote by U the set of all shift-invariant operators over K[x]. This set, endowed with the standard 
addition and composition between functions, forms an M— algebra. In the language of umbral calculus, 
U is usually named as the Umbral Algebra. 

The Pincherle derivative of an operator Q at the coordinate Xj is defined as the commutator 

Qx, = [Q,^,], (2) 

where on the right hand side of Xj acts as a multiplication operator on M.[x\, i.e. xj : p{x) t—^ Xjp{x). 

According to Theorem l2.H shift-invariant operators over polynomials are expressed in terms of derivatives. 
In particular, the multi-index derivatives and linear combinations of there are shift-invariant. 
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Regardless the last viewpoint, one looks to shift-invariant operators as convergent power series Q{x) = 
X]|^|=o"a^"' replacing the vector variable x by the gradient dx, namely = Q{dx), where 

i : M[s] — > is a mapping between the algebra of polynomials and the umbral algebra. This mapping is 
one-to-one and onto and therefore an algebra isomorphism. 

The inverse image under l : M\x\ oi a. shift-invariant operator is known as its indicator. As a further 
consequence, dxjQ{x) is the indicator of the Pincherle derivative Q'^. := Q'x-{dx) c.f. Pj. Therefore, the 
Pincherle derivative is a derivation operation in the umbral algebra hi. 



2.2 Quantum Mechanics representation of Umbral Calculus 

The main purpose of this subsection is to gather a fully rigorous description of umbral calculus in terms 
of Boson calculus associated to the second quantization. First we recall some basic facts on Quantum 
Field Theory. 

The Heisenberg-Weyl algebra !F is the free algebra generated from the vacuum vector $ by the 2n + \ 
elements a^, . . . , a~, a^, . . . , a^, id satisfying aj<i> = 0, j = 1, . . . , n and the commutation relations: 

[a7,aj:] = = [a+,a+], [aj ,a+] = Sjk'^d. (3) 

We now equip T with an euclidian inner product (-I-) such that ($|$) = 1 and the operators a^ are 
adjoint to aj: 

{a+x\y) = {x\a-y). 

The later vector space (JF, (-I-)) becomes the Bose algebra, where the elements a^ and aj (creation and 
annihilation operators, respectively) forms a basis for the n— particle space. 

The explicit expression from the basic vectors in the Fock spaces rja ^ T can be derived from the standard 
lemma of elementary Quantum Field Theory as the multi-index product of creation operators acting on 
the vacuum vector $, i.e. 

n 

and moreover, inner products between basic vector fields has the form 

In the following, we will describe the algebra of multivariate polynomials M[a;] as a realization of the Bose 
algebra !F. 

Among the infinity of possible representations of the simplest one are the multiplication Xj and the 
derivation operator dx^ ■ 

Now consider the multivariate monomials x" <E Evidently Xj and dx^ acting on a;" gives 

We now construct the umbral calculus through representations of commutation relations ([3]) by position 
and momentum operators, x'^ and Oxj , respectively such that the action of x'j, Oxj on certain polynomials 
Va{x} is analogous to the action of Xj and dx^ on the monomials. More specifically we shall search for 
x'j , Oxj and their associated polynomials Va {x) that satisfy 

X'j Va [x) = Vq+Vj [x) , Oxj Va {x) = Uj Va-^^ ix) . (5) 
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These operators can be immediately recognized as raising and lowering operators acting on the polyno- 
mials Va{x). 

In the sequel, we assume that Ox satisfy the shift-invariance property ([T]) and Ox (xj) is a nonzero 
constant. 

According to the umbral setting such an operator Ox is called a delta operator and such a polynomial 
sequence {Va : a € Nq} generated by ^ is called a Shejfer set for the multivariate delta operator 
Ox = {Oxi,Ox2,- ■■ ,OxJ; the Sheffer set {Va : a G NJ}} is called &asic if Fa (0) = and Fofe) = 1- As 
a further consequence {0'^_)^^ always exists and commute with Ox^ (c.f. [5], pp. 18). 

Using the quantum field representation (j4]) , the basic polynomials Va (x) are obtained through the action 

of fe')" :=nLiKr 

Va{x) - (6) 

which is known as the Rodrigues formula. 

A further consequence of ([5]) is the eigenproperty 

n 

Y,x'Px,Va{x) = \a\Va{x) (7) 
This fact provides a good motivation for calling 

n 

i=i 

the umbral Euler operator, i.e. an operator who measures the degree of Va] 

The properties of Sheffer-type polynomials are naturally handled within the so called umbral calculus. 
Incidentally we shall study the Sheffer operator 

^x:x"^Vaix), (9) 

where Vq,(x) is a Sheffer sequence. This linear operator is obviously invertible and its inverse ^'"^ 
corresponds to '^x^ ■ ^a{x) >— > x" ■ 

This correspondence naturally induces a correspondence between the operators Xj , dxj and Ox^ , Xj . In- 
deed, according to ([5|) the linear extension of '^x ■ x^ ^ Va{x) is an intertwining operator i.e. 

Ox, *x = ^xpx, , X'j ^x = "^xXj . (10) 

and hence it intertwines the umbral Euler operator ([5]) and the classical Euler operator Ex = xjdxj , 

i.e. 

E'^x = "fxE. (11) 

Another important consequence of (fTO|) is that '^x preserves the Heisenberg-Weyl algebra relations ([3]) . 

A formal power series is an infinite linear combination of monomials. We will use the variables x and t 
to denote formal power series whose coefficients are polynomials and 



t- X — tjXj 



to denote the formal Euclidean inner product between t and x' . 

An important consequence of the Heisenberg-Weyl representation ([3|) concerns the exponential generating 
function for Sheffer sequences Va{x), which is defined by 
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V{x,t)= ^a(x)^ = ^^te-a;')'l = expa.x')l. (12) 

|a|=0 ■ fc=0 

The above formula corresponds to a formal power series in the variables x and t. According to ()10p the 
exponential generating function (|12p may be characterized by means of V{x,t) = 5'a;(exp(t • x)). 

Let us take a close look for the characterization of Sheffer sequences on the umbral algebra. 

According to the Theorem 12.11 {^a ■ ct G Nq} is a Sheffer sequence if there is a set of multivariate 
polynomials {Oj : j = 1, . . . , n} such that 

Oj(5jT4(x) = a,V„_v, (x) (13) 

which are lowering operators. This characterization is not unique, i.e., there are a lot of Sheffer- type 
sequences satisfying equation ([13]), for a given Oj{x). We can further classify them by postulating the 
existence of the associated raising operator. A general theorem (c.f. [221 12]) states that a polynomial 
sequence {Va ■ a S Nq} satisfying the monomiality principle ([5]) with a multi-variate operator Ox = 
{Oxi , ; • • ■ I ) is uniquely determined by the polynomial vector-field 0{x) = (Oi (x), O2 (2), ■ ■ • , 0„(a;)) 
such that 

OM = 0, [dx,0,{x)]x=o^0 and [dxPkix)]x=o = for j ^ k. (14) 

The exponential generating function for Va{x) is then equal to 

V{x,t) = *3;(exp(t = exp(x • O^^t)) (15) 
where 0-\t) iO^\r),0^\t), . . . ,0-^1))- 

Furthermore their associated raising and lowering operators associated to ([5]) are given by 

Ox,^Oj{dx), (16) 
x'^=x,{0',^)-\ (17) 

We would like to emphasize that there many possibilities still for constructing operators x'j having the 
raising property. In particular, it is also interesting to consider the special choice 

Since (O^ .)~^ commutes with (c.f. [2^, pp. 18), one can easily verify the Heisenberg-Weyl commutation 
relations for the operators Ox^ and x'j defined in ([16)) and ([20| . 

One of major reasons to consider (j20p over ([TT]) is that under some special choices of O^; . , the operators 
Oxj and x'j becomes symmetric and self-adjoint on a Hilbert space c.f. [29] . 

Now we would like to point out some properties concerning the umbral operators (|16p and ([17^ : 

Since Ox is a delta operator, the operators (0^ .)~^ always exists and O'^. ., (O^ .)^^ commute with Ox 
(c.f. [a, pp. 18), i.e. 

[O;^,O.J=0=[(O;^.)-i,O,J (19) 

However the same in general does not hold between O^ ., {O'^. .)^^ and x'j defined viz formula ([TT)) (see 
Example!^. 

We would like also to emphasize that there many possibilities still for constructing operators Xj having 
the raising property. In particular, it is also interesting to consider the special choice 

■-li^AO'xX' + iO'r'x,) (20) 
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One can easily verify the Heisenberg- Weyl commutation relations ([3]) for the operators Oxj and x'^ defined 
in and pp)) . On the other hand we also get the raising property for the operators ([^0]) . However like 
for (fT7|) . the operators (^0)1 does not commute with O'^. and (0^^.)~^. 

In brief, one of major reasons to consider (PO)) instead of p?)) is that under some special choices of . , 
the operators O^j and a;^- becomes symmetric and self-adjoint on a Hilbert space c.f. [25^. 

Here there are some examples of so obtained representations of the Heisenberg- Weyl algebra: 

Example 2.1 (Forward/Backward differences) Denote by d^-' the forward/backward difference op- 
erators by 

These operators mutually commute when acting on functions, i.e. 

dt'{df,'f{x))^d^\di^f{x)). (21) 

and are interrelated with shifts T-t/jv^ 

r-H.M'f)i3L) = = [df f){x). 

Let us remark that the finite difference action d^-' acting on functions satisfies the product rule 

dh ixjfix)) = ±^ = Xjd,/{xjf{x)) + f{x + hvj). (22) 

and hence {dfj)' = tj-^jv^ . Replacing the coordinate function Xj by the operators XjT^phvj , we establish 
the duality between the finite difference operators d^'' and the "formal" coordinate functions XjT±hv. . 

We also note that the operators a:jT±/ivj mutually commute, when acting on functions 

XjT±hvjiXkT±hvtf{x)) ^ XkT±h-Vk{XjT±h-vJ{x)). (23) 

Furthermore, the commutative relations together with the duality relation \2'^) endow an alge- 

braic representation of the Heisenberg-Weyl 0^ algebra, where the "formal" coordinate functions xjTz^hvj 
represent "creation" operators dual to the "annihilation" operators 9^"' . 

Applying the Rodrigues formula 0), we obtain the following multivariate polynomials 

n 

These polynomials are the basic polynomials that satisfy the equations R is interesting to see that these 
multivariate polynomials coincide with the multi-index factorial powers {x)^^^ introduced by N. Faustino 
and U. Kahler in fESjj. 

Here we would like to emphasize that the operators XjT±hv anrf T±hv does not commute when acting on 
functions, i.e. 

T±hVjixkT±hVkfix)) = ±h6jkf{x±h{Vj + Vk)) + XkT±hVk{T±h-vJ{x)), 
T^hWj{xkT±h-Vkf{x)) = ThSjkf{x± h{-Vk - Vj)) + XkT±hv^{T±hvjf{x)) 

Another family of polynomials involving the position operator x'^ defined in i20\) can be considered. In 
particular, recursive application of the Rodrigues formula ^ yields the following family of polynomials. 

n ^ {^3'^±h^^, + ^-±h^.,XJy' l = Y[ixj±-h\ ixj±-h 

i=i i=i ^ ^ ^ 

We will come back to the Heisenberg- Weyl character of the operators 9*-' and x'j in Section when 
establishing the comparison between the Umbral version of Almansi-type decomposition and the Fischer 
Decomposition for Difference Dirac operators obtained in \25f . 
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Example 2.2 (Central difference operator) In lattice gauge theories it is common to consider the 
central finite difference operator (i.e. the average between the forward and backward differences) 

instead of the forward/backward differences d^^ . 

One among many reasons behind this choice is that contrary to Oxj are symmetric and self-adjoint 
on the Hilbert space i'2(/iZ"). 

Furthermore the Pincherle derivative for Ox^ corresponds to 

Ox, = 2 " cos\i{hdx,). 

It can be easily seen that the multivariate operator Ox — {0xi,0x2, ■ ■ ■ ,0x„) is a delta operator so 
it admits a basic polynomial sequence satisfying the equations ^ and that O'x has an inverse since 

o'x^^^^o. 

Now it arises the question how to compute (O^^)^^- Since the map l : R[ac] — > U is an an algebra 
isomorphism, the problem of finding an inverse for O'^. it is enough to find an inverse for the series 
expansion 



2k 



We however note that finding an inverse for cosh{hxj) becomes increasingly cumbersome. Nevertheless, 
there is a certain systematic in it. 

For a sake of simplicity, we impose periodic boundary conditions on our lattice such that x + Nwj = x, 
i.e. {Vivj}^ — id for a certain N E N.(i.e. a periodic lattice). In particular, when N/2 is odd, {0'x .)~^ 
corresponds to c.f. ]2S^ 

N/2-l 
fc=0 

while for N odd one obtain 

{N~l)/2 (N-l)/2-l 
k=0 k=0 

These inversion formulae will be of special interest on Section when we consider the gauged version of 
the Harmonic oscillator. 



3 Almansi-type theorems in Umbral Clifford analysis 
3.1 Umbral Clifford Analysis 

The main objective of this subsection is to introduce a weU-adapted Chfford setting. 



In what follows we retain the notation adopted in subsection 12. 21 Additionally we will denote by CZg.n the 
algebra of signature (0,n) whose generators 61,62, ... ,e„ satisfy the graded deformed anti-commuting 
relations 

{ej,ek} = -25jk, (24) 
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where {a, b} = ab + ba denotes the bracket between a and b. 

The following algebra is an algebra of radial type, i.e. 

{ [ej ,efc] ,e; } = 0, for all j, fc, Z = 1, . . . , n. (25) 

In the above definition there is indeed no a priori defined a linear space to which the generators belong. 
Nevertheless, by only using (|25|) one can already deduce many properties [23] . 



Analogously to [23] i we define Umbral Clifford analysis as the study of the operators belonging to the 
algebra of differential operators 

Alg{x;-,0,^.,e, : j = (26) 

where the operators Oxj and x'^ are defined accordingly to (fT6|) and p7p /(|20 |) . respectively. 
We define the umbral counterpart for the Dirac operator as 

n 

and by 

n 

a:' = ^e,x;. (28) 
i=i 

the umbral counterpart of the vector variable. 

The Heisenberg-Weyl character of the operators x'^ and Oxj together with the commutation relations 
([T^. and the anti-commutation rule ([M]) implies that 

n 

{x'f = -Y^{x',)\ (29) 

n 

{D'f^~Y.Olr (30) 

The operator {x'Y then corresponds to the generalization of the norm squared of a vector in the Euclidean 
space since jx'p = — (x')^l while the right hand side of ([5(7| corresponds to the umbral counterpart of 
the Laplacian, i.e. A' = — (Z?')^, and this indeed what we expected from a "true" coordinate variable 
and Dirac operator. 

Bellow there are some examples of this type of operators: 

Example 3.1 If Oxj = dxj, we recover the standard Clifford Analysis since the umbral operators D' 
and x' coincide with the standard Dirac and coordinate variable operators, respectively. Furthermore the 
continuum Hamiltonian ^ (—A + |a;p) can we rewritten as 

l(-A + |£H^l(i?2-x^), 

where the quantity ^|a;p = should be interpreted as a spherical symmetric potential operator. 

Example 3.2 If Oxj = d^'' , where d^'' are the forward differences introduced on Examvle \2.1[ the 
corresponding operator D' is given by 

D' = j2e,df^ (31) 
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The square of {D'Y' does not split the star laplacian 

71 

AVE^'^'V, (32) 
i=i 

Example 3.3 IfOxj = \ (j^h'' + '^/i"') '^''^ central differences introduced on Examvle \2.2\ the corre- 
sponding operator D' is given by 



2 



n 

D' = ^T.^^id^'+d^'). (33) 
The sguare of {D'Y split the star laplacian 



^-T.9tX, (34) 
i=i 

In Examples 13.21 and 13. 3[ the operator ~^{x')'^ describes a spherical potencial of the lattice while the 
equation — (.t')^1 = r describes a lattice sphere on the n— dimensional ambient space. 

We would like to point out that under the conditions of (|2.2D . the special choice (pil)) leads to spherical 
potentials on the lattice which possesses the symmetric property like in continuum. However the same 
does not holds under the conditions of Example 12.11 



We will come back to this problem in Section [6] when we consider the gauged version of the harmonic 
oscillator. 

3.2 Decomposition of Umbral polymonogenic functions 

Deflnition 3.1 Let D, be a domain in M" and k CzN. A function f : f2 — > Clo.n is Umbral polymono- 
genic of degree k if {D')^ f — 0. Ifk — 1, it is called Umbral monogenic function. 

Deflnition 3.2 A domain $7 C R" is a starlike domain with center if with x ^ also tx ^ holds for 
any <t < 1. 

Let r2 be a starlike domain with center 0. For any s > 0, we define the operator Is : C{fl,ClQ n) — * 
C{n,Clo,n) by 

Isfix) - / f{tx)f-^dt. (35) 
Jo 

Recalling the definition of ^>x_ in we now introduce an operator 

I's^-^Js^l\ (36) 

For any fc e N, denote 
where Ck = {-2f[k/2]\. 
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Theorem 3.1 Let fl be a starlike domain in M" with center 0. If f is a Umhral polymonogenic function 
in of degree k, then there exist unique functions /i, . . . , fk, each Umhral monogenic in 17, such that 

fix) = fi{x) + x'f2{x) + ■■■ + {x'f-\h{x). (38) 

Moreover the Umhral monogenic functions fi, . . . , fk are given hy the the following formulas: 

h = Q'^^,{D'f-'f{x) 

h-i = Q',_,{Dr-'{id-{xr-'Q'k-i{Dr-')f{x) 



(39) 



/2 - QiD' (id ~{x')^Q2iDr)---{id~{x')>^-^Qk-iiD')'^-') fix) 

/i = (id - x'Q[D'){id - {xW,{D')^) • • • (id - {xT-'QkMDT-')f{x). 



Conversely, the sum in i38\) . with fi, ■ ■ ■ , fk Umhral monogenic in Q, defines a Umhral polymonogenic 
function in fl. 

3.3 The proof of the Almansi Theorem 

Before proving Theorem 13. 11 we need some lemmas. 
Denote 

Es = sid + E. (40) 

We write E instead of Eo wlren s ~ 0. 

Lemma 3.1 (!19l ) Let x £ R" and il he a domain with ft D [0,x]. If s > and f £ C^(ri,CZo„), then 

f{x)^IsEJix)=EJJix). (41) 



Note that Lemma [3 . 1 [ remains true in the setting of Umbral Clifford analysis, since it holds componentwise. 
Since = '^^Is'i'x^ and E'^ = 'if^Es'^'^^, it follows from (gl]) that 



E'X - I'sE: = id. 



Lemma 3.2 Let fl he a domain in M". If f £ C^( fl), then 

D'E'Jix) = K+^D'fix), D'i'Jix) = l's+iD'f{x). (42) 

Proof: Recall that [O^j, = ^jfcid and [Oxj,Oxk\ = 0, i.e., 

Ox J x'l, = (5jfcid + x'^.Oxj , Ox J Oxk = Oxk Ox^ • 

By definition 

n n 



so that 



D' = }_^e,Ox,, E' = }_^x',Ox,, 

j=i k=i 

D'E' = Elk=ie,Ox,x',Ox, 

= fe=l ""ji^Jk^d + x'^^Ox-)Ox, 

~ ^jk^jOxk + X]j,fe=l ^jXf^OxjOxk 

= D' + E'D' 
= E[D'. 



12 



With the same approach, we can also show that x'D' + D'x' = —2E'^i^ and E'^x' = x'E'^^i. 
Hence D'E', ^ D'(sid + E') = sD' + E[D' = E.+iD'. Since E', = rj\ we have 



D'l: = I',+^K_,,D'I', = I',+^D'E'X = Z^+ii?'. 



As direct consequence of (UH, we find that I'^f and i?^/ are monogenic, whenever / is Umbral monogenic. 
From the definition of Q'f, in ([57)) . we thus obtain 

Q'kHi = (43) 

Lemma 3.3 Let Q be a starlike domain in M". For any monogenic function f in fl, 

{D')\x'fQ'J{x)^f{x), xen. 



Proof: From (|37|) and Lemma [3.21 we know that 

iQ',)-'=c,E'^^^^^---E'.^,E'.. (44) 

Denote g — Q'f^f- From (|43|) and (|44|) . 5 is Umbral monogenic in and 

/(x)=cfcii;'„^j^,-..i?;+ii?^5(a;). 

We need to show 

iDr{xT9{^)^CkE'^+^^yE'^+iE'.gix), (45) 

or rather 

p,)2™-i(^,2™-i)^(^) - -(m-l)!22™-ii?'„+„_i...i?;+ii?^5(2:), 
{D')^"'{x'^"')g{x) = ™!2^"£;'.+„_i---£;^+i£;^5(a;) 

for any g Umbral monogenic in il and m G N. 

We use induction to prove ((46|) . Since a;'!?' + ZJ'x' = —2E'^^2 ^-i^d E)g(x) — 0, we have 

^'(x'g(x)) = -2ii;^5(^)- (47) 
Next we show that, for any x £Vi and to G N, 



i:i'((x')2™.g(a;)) = -2m(j;')2™-i5(x); 
i?'((:r')2™-iff(:r)) = -2(a;')^("-i)i?U,n-i5(^)- 



(48) 



This can be checked by induction. Assuming that (|48p holds for to. we shall now prove it also holds for 
TO + 1. We now apply the operator x' D' + D'x' = —2E'^i^ to the function {x'Y'^ g{x): 

x'D'{{x'f"^g{x)) + D'x'{{x'f^g{x)) = -2E'^,^{[x'f"^g{x)). (49) 

By the hypothesis of induction, the first term in the left is equal to —2m{x')'^™g{x), while the left side 



equals —'^x'^"^E'^,^,^^g{x) due to the fact that 



E',x' ^x'E',^,. 



As a result, 

D'{x'^"^+ig{x)) = -x'D'{x''^"'g{x))~2E'^^^{x'^"'g{x)) 

= 2mx'^"^g{x)-2x'^"^E'^^^^^^g{x) (50) 

= -2^'""^^;/2+™5(^). 
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This proves the second equaHty of (|48|) . The first equaUty of (I48p can be proved similarly. This proves 
the identities 



Finally we apply (j48|) to prove (|46|) . It is easy to prove when m — 1. For the general case, from (|48p we 
have 

= -2m{D'Y^"'-^mx')^"''^g{x)) 

and 

= {D'fi-^-^)-^D\-2{x'fi-^-^)E',g{x)) 

= 4(m- l)(i?')''""'^"H(^')'^""''"'^^f+,„^i5W)- 

(|46p follows directly from the above induced formulas and the assumption of induction. ■ 

Now we come to the proof of our main theorem. 

Proof of Theorem [3lU Recall Hk = {f E C°°(f7,C/o.n) : {D')''f = 0}. It is sufficient to show that 

Hk=Hk-i+Tk-iHi, fceN, 
where Tk = {x')^ . Notice that Lemma 3.3 states that 

{D'fnQ'k = id. (51) 

We divide the proof into two parts. 

(i) HkD Hk-i+Tk-iH^. 

Since Hk-i C Hk, we need only to show T^_]^Hi C Hk- For any g € by ((5T|l and (|43| we have 

{D'fin-^g) = D\{D'f-'n-,Qk-M-ir'a = = 0. 

(ii) HkdHk-i+Tk-iHi. 

For any / G -fffc, we have the decomposition 

/ = (id - Tl_,Q',_,D"^-^)f + n-i{Q'k-,{D'f-'l)- 

We will show that the first summand above is in Hk~i and the item in the braces of the second 
summand is in Hi. This can be verified directly. First, 

(^')'''(id - Tfc_iQfc_i(i?')'"')/ = HdT^' - {{D' f-'n^iQ'k^i){D'f'^)f 

= {{D'f-' - {D'f-^)f = 0. 

Next, since (£>')''" V S Hi and Q^^i-ffi C Hi, we have Q'k-i{D' f^^f e as desired. 

This proves that Hk = -fffe-i + Tk-iHi. By induction, we can easily deduce that Hk = Hi + TiHi + 
■■■ + Tk-iHi. 

Next we prove that for any / G Hk the decomposition 

f^g + Tk-ifk, geHk-iJk^Hi 

is unique. In fact, for such a decomposition, applying {D')^^^ on both sides we obtain 

{D'f-^f = {D'f-^g + {D'f-^Tk-ifk 

= {D'f-'Tk-iQ'k^i{Q'k^i)-'h 

= (Q'k-ir'fk. 
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Therefore 

fk^Q'k-i{DT-'f. 

so that 

.9 = / - T,_J, = (id - n.iQ',_,{D')'^-')f. 
Thus follows by induction. 

To prove the converse, we see from (H5|) that, for any fc e N, (_D')*''+^(a;')'^i/i = 0. Replacing k by j, we 
have 

for any k > j. M 

4 The Quantum Harmonic Oscillator 

In this section, for a set of creation and annihilation operators and aj , respectively, we shall decompose 
the Hamiltonian of the classical harmonic oscillator 

into two first order operators Dj^ and Z3_ . The operators Dj^- and _D_ constitute a pair of the same 
roles as a vector operator involving position and a momentum operators, respectively. We thus can 
establish the Almansi decomposition with respect to this pair similar to our main result in the Almansi 
decomposition for the ordinary pair. 

4.1 Classical Hamiltonian splitting and Umbral Calculus 

We first consider the classical case. We set the Hamiltonian operator H = A + |a;p) and put 

As the operators satisfy the Heisenberg-Weyl commutation relations ([3]), we can split Ti. as 

H = E+- + ^id, 

where := ^t^J Hamiltonian of a field of free (non-interacting) bosons. This operator plays 

the same role of the Umbral Euler operator ([5]) . 

For D — and x = J2]=i^j^jj where ei,e2,...,e„ are the standard Clifford generators 

satisfying the anti-commuting relations 

we set D± = -^{x ^ D) ~ 12]=! ^i^J ■ Hence the operators D± and H satisfy the following relations 



{D+,D^} = -2n, [H, = -£>_, [H, D+] = D+. (52) 

Using the above relations one may easily calculate the following commutators 

[D\,Dl]=An, [n,D^_] = -2D^_, [n,D\]=2D\, (53) 

[Dl,D+\^Q, [Dl,D+]=2D^, [n,Dl]=2D+, (54) 

[D^,D\]=-2D+. (55) 



15 



From the relations ([53|) . we conclude that and Ti. are the canonical generators of the Lie algebra 

s^2(R) and hence we have a representation of harmonic analysis for the classical Harmonic oscillator (see 
e.g. [S]). On the other hand, there is such a kind of Clifford Analysis framework involving the operators 
and £>_ since 

span I iz?^, 11)^,7^1 ©span{D_,D+} 

equipped with the standard graded commutator [•, •] is a isomorphic to a Lie superalgebra of type osp{\\2) 
(see e.g. [27]). 



5 Almansi-type decomposition for the Umbral Harmonic Oscil- 
lator 

Let us consider the umbral counterpart of the Hamiltonian written in the language of Umbral Clifford 
Analysis corresponds to 

n' = \{ {D'f - {x'f ) . 

where D' and x' are Umbral Dirac and the Umbral coordinate variable defined in equations (|27p and 
respectively. Next we put 

D'^:=^{x'tD') (56) 

The representations of the Lie algebra 5/2(1^) and the Lie superalgebra osp{l\2) also fulfil in the umbral 
setting. Indeed, if we consider the Sheffer operator defined in ([9]), the umbral operators 



satisfy the Heisenberg-Weyl commutation relations ([3]). 

Moreover since D'± = '^xP±^^^ and H' = the relations (l52])-([55l) also fulfils for the operators 

_Dj_, (-D±)^ and Ti' . Now we are in conditions to formulate the Almansi-type decomposition for the 
Umbral Harmonic Oscillator. 

Here the main step consists in replace the Umbral Dirac operator D' by D'_ and the umbral coordinate 
vector variable x' by D'^. 

We take the definition of I's, Q'k as in ([9]), ((36)) and |37|) by replacing the operators x'j, E'^ and 1'^ by 
a+, (s - f )id + n' and ((s - f )id + H')'^, respectively. 

Hence we can deduce that if D'_g{x) — then 

DUiD'+r-'gix)) = -2m{D'^r-~'g{x\ D'_{{D'^f^-' g{x)) = -2(i?V)2(™-i)g(x); 
Moreover 

D'^{{D'^fg{x))=g{x). 

Let n be a domain in and € N and let (O) denote the set of all functions / : — > CZo,n such 
that {D'_)''f = 0. With the same proof as in Theorem l3.H we have the following Almansi decomposition 
in terms the eigenvectors of the Umbral Hamiltonian Ti.' . 
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Theorem 5.1 Let Q be a starlike domain in R" with center 0. If f H (f2), then there exist unique 
functions fi, . . . , fk, each in H^(fl), such that 

fix) = + D'+f2{x) + ■■■ + [D'^)''-^ fk{x). (57) 

Moreover /i, . . . , are given by the the following formulas: 

fk = Q'u-i{D'_f-'f{x) 
fk-i = Q',_,{D'_r-Hid-iD'^f-'Q',_,{DLr-^)f{x) 



f2 = Q\D'_{\d - {D'^fQ'^{D'_)^) ... (id - {D'^t-^Q',_,{D'_f-^)f{x) 

/i = (id - D'^Q',D'_){id {D'^fQ',{D'_f) ... (id - {D'^f-^Q',_,{D'_f-')f{x). 

Conversely, the sum in |57[ j, with /i, ■ ■ ■ , fk ^ ^li^)) defines a function in H^{fl). 



(58) 



6 Umbral Bridge 



We can now show that the powerful "Umbral bridge" connects continuity and discreteness. 

First we consider the continuous case. Namely, we take Ox to be the usual continuous derivative dx . 
Then the Pincherle derivative (dx-)' — [dx-,Xj] = id, D' is exactly the classical Dirac operator, so that 
the resulting Umbral Clifford algebra turns out to be the usual Clifford algebra Clo^n and hence our main 
theorem then recoveries the Almansi decomposition as in [191 . 

It is now natural to take a look for applications of the umbral calculus formalism in problems related 
to lattices. Our main objective in this section is therefore to find a Almansi-type decomposition for the 
gauged version of the Dirac Operator and for the quantum harmonic oscillator. 

First we will explain how to recover the Almansi-type decomposition in the discrete case. 



Example 6.1 (Discrete counterpart of Almansi Decomposition) Our starting point is to consider 
the difference Dirac operators introduced in J25f . In the terminology of that paper, 

i=i 

are the forward/backward versions of the Dirac operator, where d^'' introduced in Example \2.1\ and 
ei, e2, . . . , e„ being the standard basis for the Clifford Algebra Clo^n- 

By introducing the Shejfer maps 

the Almansi-type decomposition for the operators can be obtained by replacing the standard Dirac D ~ 
^j'^xj and coordinate variable operators x = ^j^j discrete counterparts 

D^ = ^^D{^^r\ and x± = *±x(*±)-i = x.t^h.,, (59) 

respectively (see \19^). 

A further consequence, we recover the Fischer Decomposition obtained in '25]. If 11^ denote the space of 
all Clifford- Valued polynomials homogeneous of degree fc, and 11 corresponds to the graded algebra 

oo 

n = 0nfc, 

k=0 
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the set 

corresponds to the spaces of Chfford-valued polynomials generated by the polynomials (s.)!^'' of degree 
|q;| = k, where = X]j=i ^j^h^'' stands the diference Euler operators introduced in [25,. 

Then the Fischer Decomposition for the Difference Dirac Operators corresponds to 

k 

= 0(^*)'-^fc-,' where X^_-=.nJnkeri?±. (60) 

We have to stress that the Fischer Decomposition obtained above does not allow to obtain a refinement 
of the discrete spherical harmonics (i.e., the null solutions of the star Laplacian/D'Alembertian (|34p ) in 
terms of null solutions of . Now we will turn our attention for the discrete case described in Example 

Example 6.2 (Refinement of discrete spherical harmonics) We shall now consider the central dif- 
ference Dirac operator acting in lattices with mesh-size h/2, i.e. 

Dh/2 ■■=\{Dl,2 + D-n„)- 

In this way, taking the central finite differences and the resulting inverse for the Pincherle derivative 
described in Example \2.S\ we build the Sheffer map 

n 

viz the symmetrized formula li20\} for the operators Xj and hence, we obtain like in Example \ 6.1\ the 
Almansi-type decomposition as well as the Fischer decomposition for the central difference Dirac operator 
Dh/2- The main advantage of the described approach we split the star Laplacian operator as the square 

of D h/2- 

Now we can construct the discrete counterpart of the spherical harmonics into the null solutions o/Z3/j/2- 
This is essentially obtained in the following way. 

Namely, if f is harmonic on the unit ballM := {x : — x^ < 1}, there exists a unique decomposition (see 
fl^ . page 1545) 

f{x) = f,{x)+xfoix) (62) 

where fi and fo are monogenic in B. Moreover, fi and /o are explicitly given by 

/i(x) = (id + xI^/^D)f{x), fo{x) = ~I^,^Df{x) 

Since the Sheffer map is an algebra isomorphism, we obtain the decomposition of discrete harmonic 
functions on the discrete ball 

by applying the Sheffer map ^ on both sides of i6S\) . Indeed we have 

^.f{x) = *./i(a;) + a:^/2%/o(x). (63) 

Moreover if we restrict the above decomposition to the discrete .sphere d{'^x^), ^^xfix) G ^^Ilfc is the 
discrete counterpart of the k— .spherical harmonic polynomial f{x) while 'i'xfiix) G "^x^k is the discrete 
counterpart of the k— spherical monogenic polynomial fi{x). 
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From the definition, we can also see that contrary to the operators , the central difference Dirac 
operator Dh/2 is not a 'local operator' since the shift operators Tj^h^. when acting on some lattice 
functions, they not only concern with the nearest neighbor points but also all the points contained in the 
direction h^^j. So we can call the combination I?;i/2 a 'quasi- local' operator. The quasi- local is the price 
for getting the above mentioned nice description of discrete spherical monogenics. 

The next important example concerns the Gauged version of Harmonic Oscillator written in terms of the 
operators described in the above example. 



Example 6.3 (The gauged version of the Quantum Harmonic Oscillator) Now we consider the 
Gauged version of the Quantum Harmonic Oscillator. 



This local operator is defined on the uniform lattice with mesh-width h and corresponds to a splitting in 
terms of 'quasi-local' operators. 

In order to describe the eigenf unctions of the gauged Hamiltonian, we turn now our attention to the 
continuous harmonic oscillator described in subsection\4-.l\ 



In the continuum, a straightforward computation shows that 

Va{x) = exp(-\xf)Ha(x), 

are eigenf unctions of TL with corresponding eigenvalue X — |a| -t- ^, where IIa{x) stands the classical 
multivariate Hermite polynomials. 

For these eigenf unctions, we also obtain following raising and lowering properties for the set of operators 

Va {x) = Va+vj {x) , sq {x) ^ aj Va-vj (x) . (64) 

Moreover, the exponential generating function for Va (x) is given by 

Vix, t)^Yl = exp [-1 {\x\' + lip) +V2x-t 

The corresponding eigenf unctions and exponential generating function for the gauged Hamiltonian Tih 
are then obtained viz the Sheffer map \61]) . These are closely related with the Kravchuk polynomials of 
discrete variable introduced in JSl}/ . 

Furthermore for the spaces of Clifford valued polynomials generated by "^xVaix) we obtain the following 
version of the Fischer Decomposition 

k 

@{xu/2y^Jl+-,{^), (65) 
where ^^H+_^in) := {/ G H : Hhf = (fc - j + §) /} n ker (A./2)-- 



In the last example, the Fischer decomposition results as a direct consequence of the Howe dual pair 
technique (see pSj). 

Moreover, the following theorem appears as a generalization of (|65p. allow us to split the eigenvectors of 
the Gauge Hamiltonian Tih in terms of null solutions for the operator {Dj^/2)-- 
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Theorem 6.1 (Almansi-type decomposition for the Gauge Harmonic Oscillator) Letfl be a star- 
like domain in M" with center 0. // / satisfies 

(Av2)-/(a;) = 0, and Hhf = [k + ^)f{x) (66) 

in n (/i/2)Z", then there exist unique functions /i, . . . , fk G ker(Z?^y2)-, such that 

fix) = hix) + [Dy,,2)+f2{x) + ■■■ + {D,-,,.2)X-^fk{x). (67) 

Moreover the discrete monogenic functions fi, ■ ■ ■ , fk o,re given by the the following formulas: 



fk = Q'k-iiDh/2)--'fix) 



fk-1 = QUPV2)- (id-(^V2)r Qfe-i(^V2)^')/( 



f2 = Q'i(A./2)-(id-(i?V2)^g2(A./2)-)---(id-(i^V2)+ Qfc-l(A^/2)-"')/(:r) 

/l = (id - (A^/2) + Qi(A^/2)-)(ici - ( A./2 )^ Q2 ^2 )- ) ' ' ' (id - (i^V2)+"'Qfe-l ^.72)-" V(a:) • 

Conversely, the sum in |6'7| j, the functions /i, . . . , /fc satisfy the equations id 6]) in flh/2- 
With the above resuhs show the powerful Umbral bridge unifies continuous and discrete. 



7 Concluding Remarks: 



Umbral Calculus is an analysis of correspondence between continuity and discreteness in the commutative 
field related to polynomials. In the article, we bring Umbral Calculus into the non-commutative field 
through the Umbral Clifford analysis and we further consider functions far beyond polynomials. We 
study the discrete version in the Umbral Clifford analysis by keeping the symmetry of phenomena with 
the price of non-commutative. The Almansi decomposition is established in Umbral Clifford analysis, 
which concerns the theory of the powers of Laplacian. In particular, it gives the Fisher decomposition of 
polynomials, which is a start pointing for the spherical harmonics and special function theory as a whole. 
One can further study the theory of Fourier transformations and wavelets. Recently the theory of power 
operators of Laplacian has found their applications in the conformal Geometry in pure mathematics (see 

m)- 

As a framework of the unification of continuity and discreteness, the strong power of this framework lies at 
its constructive isomorphism. This makes it possible to construct the associated equations and solutions 
in the discrete setting starting from the equations and their solutions in the continuous setting in Physics. 
One may further bring this technique to the field of Geometry by studying the fruitful equations in the 
continuous setting through the associated discrete version. 
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